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Oxidation reaction of CO on a single platinum crystal is a reaction-diffusion 
system that may exhibit bistable, excitable, and oscillatory behavior. We stud- 
ied the effect of a stochastic signal artificially introduced into the system through 
the partial pressure of CO. First, the external signal is employed as a turbu- 
■ lence suppression tool, and second, it modifies the boundaries in the bistable 

transition between the CO and oxygen covered phases. Experiments using pho- 
' toemission electron microscopy (PEEM) together with numerical simulations 

performed with the Krischer-Eiswirth-Ertl (KEE) model are presented. 



1 Introduction 

^SJ ' Chemical reactions in extended media are one of the characteristic systems 

■ where pattern formation outside equilibrium can be observed. The nonlinearity 

of such reactions and the diffusion of its components induce the formation of 
complex patterns. The Belousov-Zhabotinsky reaction and the catalytic oxi- 



CN ■ dation of CO on single crystal platinum surfaces are chemical examples where 



spiral waves, target patterns, or chemical turbulence have been studied [T]. Such 
systems have been extensively employed because it is relatively easy to control 
the chemical and physical conditions which are responsible for the formation of 
the spatiotemporal structures [2]. 

Natural systems are unavoidably subject to random fluctuations (noise), 
?H ' stemming from either environmental variability or thermal effects. These unde- 

sirable perturbations mask and affect the deterministic dynamics of the system. 
Such fluctuations gain importance far from equilibrium [21 H] because they can 
interact with the nonlinearities of the system, giving rise to a rich variety of 
behaviours [5] . The effects of noise in extended active media [6] are particularly 
interesting. Chemical reactions are one of the most convenient systems used 
to study this type of phenomenology. Both Belousov-Zhabotinsky and chlorine 
dioxide-iodine-malonic acid (CDIMA) reactions are good examples where the 
effects of noise have been studied extensively. The photosensitive version of 
such reactions is particularly suitable for the introduction of controlled external 
fluctuations through a stochastic illumination [3 [HI [S] j where the characteristic 
time and length can be easily tuned by a computer. The most significant re- 



1 



suits have been obtained for noises with spatial structure. However, interesting 
results have also been observed with homogeneous noisy illumination [10] or 
applying an stochastic electric field [TT] . 

Heterogeneous catalytic CO oxidation is well known for the wealth of vari- 
ety of pattern formations that can develop on the surface of a single crystal. 
Depending on the partial pressure of the species inside the reactor chamber, 
different types of patterns have been observed: spiral waves, target patterns, 
standing waves, and turbulence [H]. This reaction is not photosensitive, but 
the pressure of the reactants can be externally controlled |13| : in particular, a 
stochastic signal can be artificially introduced. There is no spatial structure in 
the external signal and therefore the noise is homogeneous. Previous experimen- 
tal and analytical studies in bistable reactions show how noise can anticipate 
the transitions between both bistable states in CO oxidation in Ir(lll) [Ml US], 
and how local fluctuations become important |16j . Theoretical predictions of 
stochastic resonance in the catalytic oxidation of CO in Pt(llO) have also been 
reported [T7j . 

The effects of external noise can be interpreted as a control mechanism of the 
different types of patterns, and by tuning the characteristics of these external 
fluctuations, transitions and new states can be achieved. The control of chemical 
turbulence is a problem tackled from a number of different angles [U [TB] . In the 
CO-Pt system for instance, attempts have been made with feedback techniques 
[TOl [2U], local control [31], periodical global forcing ^22], and resonant global 
forcing [T3] . 

Two different scenarios of noise effects are considered here: first, we study 
the elimination of turbulence through stochastic forcing, and second, we focus 
on the noise related modification of the boundaries in the bistable transition 
between the CO and oxygen covered phases. Both experimental studies are 
contrasted with the results of numerical simulations performed with the KEE 
model for the CO oxidation in Pt(llO) ^Sl- Finally, we discuss the results and 
compare them with previous control strategies for the catalytic CO oxidation. 

2 Experimental setup 

To perform the experiments, a single crystal Pt(llO) is mounted inside a con- 
stant flux ultra high vacuum (UHV) chamber, where CO and oxygen are con- 
stantly added into, while the total pressure is kept constant by continuously 
pumping the entire system. A PEEM [24] allows us to spatially resolve the lo- 
cal work function of the adsorbate species and to record it with a CCD camera 
at a rate of 25 full frames per second. A meaningful modulation of the CO 
partial pressure was introduced into the system as a computer controlled source 
of random fluctuations following: 

Pco[U)^P%o{^ + ^,-^N) (1) 

where ipi is a random number between -1 and +1 that is regenerated by the 
computer in time intervals of At (ti — ti-i = At). AiV is the percentage of 
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Figure 1: (Color online) Histograms of the pco{t)- a) Generated signal, b) 
signal in the dosing system, and c) mass spectrometer signal inside the reaction 
chamber (Gaussian distribution overplotted) . 

CO partial pressure that is added to the base value. The combination of the 
parameters At and AA'^ determines how fast and how much pco{t) changes, or 
in other words, the magnitude of the multiplicative noise present in the reac- 
tion. The two parameters are then used by the computer to control the CO gas 
dosing valve. Due to the solenoid valve effective reaction time, and the low pass 
filtering between the dosing system and the chamber, the generated rectangular 
signal is integrated by the system, resulting in a Gaussian-like distribution in 
the solenoid valve and in the reaction chamber, as shown in Fig. [TJ Hence, 
the CO partial pressure in the reaction environment has a fluctuating value 
between p^Q + AN and Pqq — AN, homogeneous in space, characterized by 
a histogram like the one depicted in Fig. [l]b which resembles a Gaussian noise. 



3 Reaction model 

In this work we use the KEE model [23], which has been previously employed 
to study CO oxidation in surfaces of Pt(llO) [25l[26]. The model consists of 
three coupled equations for the CO and oxygen coverage, respectively u and 
and the fraction of the surface found in the nonreconstructed Pt 1 x 1 phase, 
corresponding to the third variable w. The equations read: 

D\7^u + kiScoPcoC^ - u^) - - k^uv, (2) 
k4P02[so.ixiw + So, 1x2(1 - w)]{l - u-vf ~ kjuv, (3) 

'■^( i.^pW) -")- « 

For a brief explanation of the model and the parameter values used here, see 



dtu = 
dtv = 

dtW = 
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Table I. Such parameters correspond to previous numerical studies of the same 
experimental system under similar conditions [25j . The values of the partial 
pressures of CO and oxygen are respectively pco and po^ , and they are specified 
for each case studied here. Several variants of this model have been previously 
employed to study different types of bifurcations observed in the CO oxidation 
or the dynamics of spiral waves [57j. The effect of the temperature can also 
be considered, assuming an Arrhenius-type dependence on the temperature of 
the model parameters |26) . Here, we keep the parameter corresponding to the 
temperature constant T = 542.35 K. 
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Figure 2: Phase diagram of the KEE model keeping the values of table I con- 
stant. This diagram does not take into account the diffusion. Regions of oscil- 
lations and bistability are marked respectively as black and dark gray areas. 

The phase diagram of the KEE model for both partial pressures {pco and 
^02)1 keeping the rest of the parameters constant, is shown in Fig. [2] It cor- 
responds to the system without diffusion, thus there is no information about 
spatial patterns. As expected for large values of the pressure of oxygen and low 
pressure of CO, the oxygen covered phase is stable. On the other hand, for large 
values of the pressure of CO and low pressure of oxygen, the CO covered phase 
is also stable. For intermediate values of both partial pressures (typically with 
the pressure of oxygen three times larger than the pressure of CO), both phases 
compete. While for large pressures this competition gives rise to bistability and 
hysteresis, for intermediate pressures it also induces oscillations. In this region 
of periodic alternation, not only global oscillations but turbulence are observed 
depending on the values of the parameters. For low pressures of both CO and 
oxygen, a continuous transition between both phases is observed (dashed line in 
Fig. [21 This diagram qualitatively agrees with the experiments of CO oxidation 
on Pt(llO). 

In order to mimic the experimental conditions of random external modula- 
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tions of the partial pressure of CO, we add a stochastic variable to the value of 
the pressure: 

Pco^p"co{^ + Ut)) (5) 

where the variable ^e(i) is a Gaussian spatially homogeneous white noise with 
null mean < ^e{t) >= and variance < ^e{t)£,e{t') >= 2a^S{t — t'), where cr^ is 
the intensity of the noise. 

It is known that the surface of the crystal is not completely homogeneous 
and there are heterogeneities which affect the propagation of the surface waves 
qualitatively. Such heterogeneities are consequence of several experimental fac- 
tors (facetting, subsurface species or cristal surface defects) and their activity 
may depend on several processes like the local concentration of CO and O in an 
unknown way. In order to consider these internal fluctuations, we introduce an 
additive noise term £^i{x,t) in the first equation of the model Eq.([2]). This vari- 
able is Gaussian white noise, distributed in space with null mean < ^i{x, t) >= 
and variance < £,i{x,t)^i{x' ,t') >= 2ei6{x,x')S{t — t'), where Ci is the intensity 
of the noise. We note that these internal fluctuations are not related with the 
external experimental noise which has been already introduced above as mul- 
tiplicative noise. Since we are not going to do a systematic study of these 
internal fluctuations, we are keeping the intensity of this noise low and constant 
{ei — 0.0005) for all of the numerical simulations following. They will act as 
precursors in the simulations during transitions. 

In summary, we introduce two types of random processes: first, a homoge- 
neous random noise through the parameter pco , mimicking the external exper- 
imental noise introduced by the CO pressure, which is the main object of study 
in this manuscript. Second, a spatiotemporal noise, that takes into account 
different experimental sources of randomness. 

We integrate Eqs. mousing a finite differences scheme and the Heun method 
|28| for stochastic partial differential equations. Two-dimensional numerical 
simulations are completed using a finite differences scheme with a spatial and 
temporal discretization of dx — 2/im and St = 0.001s respectively. We employ 
periodic boundary conditions in all of the simulations. 

4 Turbulence suppressed by noise 

Spatiotemporal chemical turbulence spontaneously develops on a platinum crys- 
tal surface for a window of values of the reaction parameters (reactants, partial 
pressures, and temperature). Such turbulence is characterized by the contin- 
uous creation and annihilation of phase defects on the surface [29j. Although 
it corresponds to a complex state, a characteristic frequency can still be ob- 
tained. Such turbulence can be effectively suppressed by a noise modulation of 
the reactants, namely pco{i)- 

Fig. [3] represents the corresponding spatiotemporal patterns seen when dif- 
ferent noise intensities are applied to a turbulent regime. In the absence of 
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Table 1: Parameters of the reaction model. 



Parameter Value Meaning 



fci 3.14 X lO^mbar ^ Impingement rate of CO 

k2 10.21s~^ CO desorption rate 

ks 283.8s- 1 CO Reaction rate 

ki 5.86 X 10^s^^mbar~^ Impingement rate of O2 

ks 1.61s~^ Phase transition rate 

SCO 1-0 CO sticking coefficient 

sOixi 0-6 Oxygen sticking coefficient on the 1x1 phase 

S0ix2 0-4 Oxygen sticking coefficient on the 1x2 phase 

Mo, du 0.35, 0.05 Parameters for the structural phase transition 

D 40/im^s~^ CO diffusion coefficient 

Uref 0.3358 Reference CO coverage 



external noise, spiral turbulence develops freely (see 2D images in Fig.([3la)). 
For large intensities of the fluctuations (above 0.2 % Pco), the noise completely 
annihilates the reaction dynamics; the whole system then follows the noisy sig- 
nal Fig. [SJd. For intermediate noise intensity, the system is driven to a partially 
suppressed state featuring an alternation between asynchronized and homoge- 
neous global oscillations Fig. [3]i. It corresponds to a competition between 
the turbulence behavior of the surface and the global synchronization produced 
by the noise. This alternation was observed to be regular and periodic for a 
particular range of parameters Fig. [3}:. Remarkably, the frequency of such re- 
sponses are different than the characteristic frequency of the spiral turbulence. 
Finally, noise reaches values that completely annihilate the reaction dynamics; 
the system then follows the noisy signal. 

The natural frequency of the system is affected by the noise. Fig. [3] shows 
how the periodicity of the local oscillations changes with A7V visible within the 
snapshots and the space-time plots, and compares the initial spiral turbulence 
of region I with the periodic induced patterns within region II, with the noise 
controlled region III for a particular fixed value of At. Region III has no local 
periodicity whatsoever and reflexes the imposed random oscillations. 

Region II of Fig. [3] encompasses the range of noise parameters in which 
the original spiral turbulence is not completely suppressed. Interaction with 
noise drives the active media through transient states, alternating between spa- 
tiotemporal chaos and a lack of patterns. Furthermore, for certain parameters 
these cycles may become regular. At the bottom left of Fig. |3] suppressed pat- 
terns alternating with re-activated turbulence in regular cycles of 4 seconds can 
be seen. The noise parameters corresponding to this induced behavior draw 
a narrow stripe in the turbulence suppression map (seen red bar in Fig. [3]). 
Along that stripe, the interplay between excitability, refractory time, noise, and 
diffusion results in time regularity. 
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Figure 3: (Color online) PEEM snapshots A of the Platinum surface and space- 
time plots B, for different values of A7V applied on an initial turbulence state 
of the chemical reaction. The spacetime plots represent the reaction pattern 
development in time along the white bar depicted on A. The corresponding pixel 
image intensities are plotted on C. Reaction parameters are T — 540if , = 
1.3-10~^m&ar, = 6.6-10~^m6ar', and At =0.5 s. Three regime differentiated 
areas are pointed out (I, II, and III), showing turbulence (a), noise-controlled 
oscillations (b), periodic alternation (c), and alternating patterns (d), which are 
further commented on the text. 



Larger or faster modulations than those showed in Fig. [3] are not considered 
here because of the limitations of the experimental setup (experiments showed 
strong disagreement between the computer values and the real state in the 
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chamber for values At < 0.1 or AN > 3 % Pco)- Therefore this study project 
was specifically focused on the smallest fluctuations with a relevant impact. 

PEEM KEE 




Figure 4: Compared space-time plots corresponding to the turbulence sup- 
pression map regions I, II, and III; experiments (PEEM) vs. simulations 
(KEE). Simulations are performed with the values of the partial pressures: 
PO2 — 1.3 • 10^"' mbar and = 4-10^^ nibar in a system of size 0.2 x 0.2 mm^. 
From top to botton the corresponding noise intensities are a= 0.0 % Pco^ 
0.25 % p"co and a= 2.5 % p^Q. 

The series of space-time plots depicted in Fig. |4] show the typical surface 
dynamics of the three regions (I, II, and III) adressed in Fig. [3] compared with 
numerical simulations; one can see that the plots based on the KEE model qual- 
itatively agree with the experimental results. Wave turbulence (I), alternation 
of global oscillations and wave turbulence (II), and dominance of the noise (III) 
are observed in both cases. 

The results of the numerical simulations agree and confirm the experimental 
results. Under parameter values corresponding to turbulence in the KEE model, 
the introduction of a global stochastic signal may induce the global coordinating 
of the whole system. For intermediate values of the external signal a competition 
of the global coupling and the deterministic turbulent dynamics is observed, 
giving rise to intermittent turbulence and a number of intermediate patterns. 
For some noise intensities, periodic oscillations are observed and even standing 
waves or pacemakers are temporarily created by the global noise. 

In the numerical simulations, the frequency of the spatiotemporal patterns 
also decreases with the noise intensity. The patterns of the chemical turbulence 
can be characterized by large frequencies. The global noise not only couples the 
whole system but also reduces the frequency of the patterns. 
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5 Transitions anticipated by noise 



Further experiments were focused on the bistable transition from a mainly oxy- 
gen covered Pt crystal to a mainly CO covered state. These experiments start 
with a certain amount of oxygen in the reaction chamber (j>02 = 1-1 • 10^^ mbar) 
with no CO. Once the surface is oxygen covered, the is kept constant and 
then the pco is increased while adding additional fluctuations to it. For a large 
enough CO partial pressure, all of the adsorbed oxygen reacts away releasing 
CO2, the sample is completely CO covered and the transition is complete. This 
process is clearly discernible in the PEEM images (top of Fig. O: it is seen 
as the change from a uniform dark area to a uniformly bright one through a 
number of growing bubbles. Due to the introduction of the additive noise in 
Eq. [21 the KEE model casts similar transition patterns (bottom of Fig. [5]). 
The intensity of such additive noise has been tuned to qualitatively reproduce 
the experimental patterns shown in Fig. O and kept constant for all of the 
simulations. Without such noise the system basically behaves homogeneously 
and no bubbles would be observed. 

The same steps are followed backwards for the back-transition to the initial 
CO covered state. Nevertheless, the critical pco is lower in this case, which is 
explained by the asymmetric inhibition of the reaction: the oxygen covered Pt 
does not behave in the same way as the poisoned CO covered surface [30]. There 
is an inherent hysteresis loop in this bistable transition due to the asymmetry 
of the reaction. Fig. [SJ'V describes how the mentioned hysteresis loop shrinks as 
the noise applied is larger. 

Equivalent numerical simulations in the bistable region reproduce the main 
features of the hysteresis loop observed in the experiments. Starting from an 
oxygen covered phase, the pressure of CO is slowly increased inducing the tran- 
sition to the CO covered phase. The posterior reduction of the CO pressure 
shows a lower value for the transition, giving rise to the hysteresis loop. Two 
numerical realizations with different noise intensities of this loop are shown to- 
gether in Fig. [6j3. Larger noise intensities imply smaller hysteresis loops. This 
reduction of the size of the loop is plotted in Fig. [6p. 

6 Discussion 

When the CO-oxidation on Pt system shows spiral turbulence and global noise is 
introduced, this noise behaves as an overall coordinating factor that introduces 
a certain spatial coherence to otherwise spatially discoordinated oscillations. 
In this way, chemical turbulence can be suppressed. The cases of intermittent 
turbulence suggest that the global coordination is not an absolute one for in- 
termediate noise intensities, and that the time needed for the active medium to 
recover plays a role in the alternation periodicity. 

In the numerical simulations performed adding external noise to the KEE 
model, the effects of the fluctuations are qualitatively equivalent to the ex- 
perimental case. Low noise does not modify the dynamics of the system and 
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Figure 5: PEEM snapshots (top) compared with numerically generated images 
(bottom). Different stages of the bistable transition, oxygen covered to CO 
covered can be seen. Experimental parameters are: T = 509 K, po^ = 1-1 • 10""* 
mbar, At= 0.25 s, and AA^ = 0.5 % Ppo- Parameters used for the simulation 
are: po^ = 4.5 • 10"'*, pco — 9.08 • lO"'"^ in a system of size 0.4 x 0.4 mm'^, and 
internal noise Ci = 0.0005. 

large noise eliminates turbulence for similat noise intensities. For intermediate 
noise intensities, alternance between global oscillations and complex dynamics 
are also observed mimicking the experimental results. The accuracy of the nu- 
merical simulations allows the observation of a richer phenomenology for these 
intermediate intensities. Moreover, other interesting patterns induced by the 
noise are also witnessed in the simulations, like the formation of clusters or 
the generation of local pacemakers during several oscillations. Two examples 
of such small pacemakers in the same numerical simulation can be seen in Fig. 
[71 where the onset of the turbulence is delayed by the noise (too weak to com- 
pletely supress the turbulence) and two transitory pacemakers appear instead 
of the turbulence. Similar types of patterns have been previously observed in 
the same system under global feedback 1 23] . 

The control of chemical turbulence has been previously studied with other 
techniques implying global resonant periodic forcing [13], feedback mechanisms 
[19], or local reaction modifications via microstructures or laser [21]. The forc- 
ing amplitude, percentage of CO partial pressure, needed to effectively suppress 
turbulent states was reported to be 10-20 % for the periodically forced reaction 
|22| . Closer analysis of the initial turbulent dynamics, determining the natural 
frequency of the system, and performing periodical global forcing at a resonant 
frequency [13] allows us to reduce that needed percentage in two orders of mag- 
nitude. As shown here, the same order of magnitude (0.1-0.2 %) is enough for 
fast random modulations to reach total entrainment. It is remarkable that the 
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Figure 6: Panel a) shows the critical pco of both transitions, from oxygen 
covered to CO covered (solid line) and from CO covered to oxygen covered 
platinum (dashed line), for a value of applied noise. Reaction parameters are: 
T = 509 K, P02 = 1-1 ■ lO""* mbar, and At = 0.1 s. Panel b) shows the 
hysteresis loop for two noise intensities cr= 0.09 % pco and a— 0.28 % Pco- 
Panel c) represents the corresponding numerical simulation plotted in terms of 
the simulation parameter a. Numerical simulations have been completed with 
PO2 = 3.5 • 10""* mbar in a system of 0.2 x 0.2 mm^. 

latter is independent from our previous knowledge of the system. 

The intensities of the external noise which eliminates the turbulence in the 
numerical simulations are also quite low (0.2-0.5 %) and they are comparable 
to the case of resonant forcing [31 . 

The fluctuations also have an influence on the transition between two stable 
states, anticipating a bistable transition of the system and shortening the intrin- 
sic hysteresis loop. Using a mechanistic heuristic model, the bistable transition 
anticipation is conceptually equivalent to the motion of a particle localized in 
the higher minimum in a tilted double- well potential [32] . in which external 
noise adds energy to the system to save the potential barrier. For moderate 
noise it has been shown that the reaction mechanism of the system more easily 
jumps the barrier that separates the two steady regimes. For large noise, this 
barrier would become so low in front of the fluctuations that the system would 
always be at one of the two local minima of the potential or transiting between 
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Figure 7: Space-time plot of the evolution of the system under turbulent condi- 
tions and with external stochastic forcing. Snapshots of the system are plotted 
in the pannels. Panel A shows the initial global oscillation typically observed 
in the numerical simulation. Panel B shows the growth of a single spot which 
finally disappears. Panel C shows the generation of a pacemaker due to the 
intermediate amplitude of the external global noise. Panel D shows the tur- 
bulent patterns. Values of the partial presures are: po-^ = 1-3 • 10^^ and 
Pco = 4.81 • 10~^. Internal and external noises are respectively e.^ — 0.0005 
and a = 0.3%pco- The size of the system is 0.2 x 0.2 mm?. 

them. 

Defect-induced fronts were experimentally observed during the anticipated 
bistable transitions. When the transitions were completed and the crystal sur- 
face showed a full CO coverage, deep drops of the pco signal may led to the 
emergence of spirals and target patterns at some defects of the surface. These 
oxygen structures propagated along the surface for a short time before they 
disappeared. 

Besides the spatiotemporal noise-induced patterns, experiments performed 
with noise also featured an increase in the resulting amount of CO2 produced in 
the CO oxidation. A mass spectrometer measured the CO2 present in the UHV 
chamber while chemical turbulence developed, with and without modulating 
the Pco- For the same amount of reactants involved in the reaction, integrated 
in time, the CO2 production rate was higher with a fluctuating CO partial 
pressure. Higher productivity thanks to fluctuating parameters is a fact already 
observed in other environments, namely the ozone production in atmospheres 
of urban areas [32]. This fact calls attention to the possible catalytic effect of 
noise, which might have some impact in the increase of chemical productivity 
for some other systems. 

In the context of such real-world systems like marine ecosystems [331 131] 
or ozone production [32], in which stochasticity processes combine with peri- 
odic annual or daily cycles, it seems worthwhile for the future to combine the 
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two global modifications: periodic and noisy modulations of the CO oxidation 
reaction parameters. 

In conclusion, we have shown that a stochastic external global signal in- 
troduced in the CO oxidation on Pt(llO) through the partial pressure of the 
CO can eliminate spatiotemporal turbulence with low characteristic amplitudes 
of the noise, comparable and even smaller than in equivalent studies with the 
amplitudes of periodic forcing. Furthermore, this artificially noisy signal can 
anticipate phase transitions in bistable conditions by reducing the reaction in- 
trinsic hysteresis loop. 
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